arXiv:1509.02288vl [hep-lat] 8 Sep 2015 



ELSEVIER 


Available online at www.sclencedlrect.conn 

*•#' ScienceDirect 

Nuclear and Particle Physics Proceedings 00 (2015) 1-0 


Nuclear 

and 

Particle 

Physics 

Proceed 

ings 


ABOUT MONOPOLES IN QCD. 


Adriano Di Giacomo 

Dipartimento di Fisica, Universita’ di PISA, 3 Largo B. Pontecorvo, 56127- Pisa, Italy. 


Abstract 

The hypothesis is analysed that the monopoles condensing in QCD vacuum to make it a dual superconductor are classical solutions 
of the equations of motion. 
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Introduction. - 

Confinement of colour is an experimental fact. No 
quark has ever been observed neither in ordinary matter 
nor as a product of high energy particle reactions. 

An upper limit has been put by experiment to the ratio 

of the abundance in nature of quarks to the abundance 

of protons — < 10“^^ , 15 orders of magnitude smaller 
than the expectation — 10“'^ estimated in the Stan- 

dard Cosmological Model on the hypothesis that quarks 
propagate as other particles (ll]. 

Similarly the inclusive cross section erg for producing 
a quark in a high energy reaction initiated by hadrons, 
which should be of the same order of magnitude as the 
total cross section o-j, has an upper limit ^ < 10“'^ . 
The natural explanation of an inhibition by a factor 
10 is that confinement is an absolute property, i.e. 
that rig and o-g are strictly zero, protected by some sym¬ 
metry. Deconfinement is then a change of symmetry. 
This symmetry can not be a subgroup of the gauge 
group, since in the Hilbert space with quarks the rep¬ 
resentation is faithful and gauge invariance is believed 
to be valid both below and above the transition. It is 
most probably related to topological excitations on the 
spatial boundary at infinity, as happens in well known 
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models of statistical mechanics ind + \ dimensions: in 

the 1 - 1-1 dimensional Ising model the excitations are 
kinks ||2l, in the 2 - 1-1 dimensional X-Y model the exci¬ 
tations are vortices |13l, in the 3 H-1 dimensional compact 
t/(l) gauge model on the lattice they are monopoles |01 
ll^ . The nature of the dual excitations is topological and 
is determined by d: they are monopoles for d = 3. 

The mechanism by which QCD confines color is not 
easy to study, since there is no precise understanding 
of QCD at large distances. An attractive idea is that 
vacuum is a dual superconductor, due to condensation 
of magnetic charges 0] 10]: dual means that electric 
and magnetic are interchanged with respect to ordinary 
superconductors. The chromo-electric field acting be¬ 
tween aq - q pair is constrained into an Abrikosov flux 
tube, producing a potential proportional to the distance, 
in the same way as is the magnetic flux in an ordinary 
superconductor. Such a mechanism includes both the 
idea that deconfinement is a change of symmetry, and 
that the degrees of freedom involved are topological ex¬ 
citations at the spatial boundary. 

The mechanism has been tested in detail in A = 2 su¬ 
persymmetric QCD and in compact t/(l) gauge the¬ 
ory on the lattice|0]0. Unfortunately real QCD is not 
supersymmetric, and the compact t/(l) gauge theory is 
probably not even a field theory. When it comes to phys¬ 
ical world things become more complicated. 

Lattice is a non perturbative formulation of QCD allow¬ 
ing to explore large distances. Monopoles can be de- 
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tected in the lattice configurations by the recipe of Ref. 

Having in mind the gauge group S U(2) for the sake 
of simplicity, the idea is to fix a gauge [Abelian Pro¬ 
jection] and to look at the flux through the surface of 
elementary cubes of the diagonal component (3 compo¬ 
nent in colour space) of the magnetic field. A monopole 
is detected as a Dirac string, i.e. as a 2n excess of 
the flux through a face of the cube ifioll . Both the 
positions and the number of the observed monopoles 
are gauge dependent. A particular gauge, the Maximal 
Abelian Gauge |@], has been privileged in the literature 
because the field configurations are dominated by the 
monopoles ifTlll . However an attempt to determine from 
the monopoles observed in lattice configurations an ef¬ 
fective action with a minimum at a non-zero value of the 
monopole field was not successful lIT^ ifoll . An alterna¬ 
tive approach was to construct a magnetically charged 
operator and use its vacuum expectation value as an or¬ 
der parameter for condensation lEl iS. The result was 
positive, showing that in the confined phase vacuum is a 
dual superconductor and normal in deconfined, but the 
result was a kind of average over the Abelian Projec¬ 
tions. 

In this work we have a different approach: we look for 
monopoles as classical stable field configurations and 
analyse their observable effects, and their possible con¬ 
densation. 


2. Monopoles in Gauge Theories. 


Monopoles are topological configurations in 3d, i.e. 
non-trivial mappings of the sphere 5 2 at spatial infinity 
onto S0{3)IU(l). They are stable against local fluctua¬ 
tions. They can coexist if their topologies are compati¬ 
ble El. The prototype monopole was discovered in the 
SO(3) Higgs model El El . 






fp- = ^ . For a static solution do = 0. In particular 
doAo - 0 , which we shall assume as a choice of gauge. 
The static equations of motion are 

^ ji ~ ^ ^ OiAo 

DiDiXo = (1) 

- yuE = 0 


In presence of a non zero Higgs field a solution exists 
with Aq = 0 obeying the equations El El 


Djdji = A Di6 

DiDi$ + - fp) = 0 


In QCD there is no Higgs field, <5 = 0 and Eq’s ([TJ read 


- ~8‘^o A DiAo 

DiDiXo = 0 


(3) 


The minus sign in Eq([3]l comes from Minkowski met¬ 
rics. In terms of A 4 - iAo Eq’sQ have the same form 
as Eq’s (|2]i with <1> replaced by A 4 |18| and d = 0. 
Monopoles in QCD have the same form as those of Ref. 
1 1611 EtIi . with A 4 playing the role of the Higgs field. 
A - 0 implies that P - ±iJ [Self-duality or anti-delf- 
duality][ 19] E is the euclidean electric field, i.e. i 
times the physical electric field. Explicitly the solutions 
are, in the ’’hedgehog” gauge : 


A'l ^ eain-[l - m)] 
gr 

A" = 

8 


(4) 




_±_ 

sinhf 


//(^) = coth^-- 


(5) 


Here ^ = pr; the upper index runs on colour, the lower 
index on space-time components. The solution is a 
mapping of the surface at infinity S 2 onto S0(3)/U{l). 
K(0) = 1, H(0) = 0, K(oc) = 0, H{oo) = 1. There is 
one free parameter, ju, which is related to the value of 
A 4 , the effective Higgs field, at r = 00 . 

The magnetic field B‘‘ at large distance has the form 


£> ■ 



( 6 ) 


The direction of A 4 in color space can be made r in¬ 
dependent by a time independent gauge transformation 
U{r) transforming to the so-called Unitary Gauge. The 
transformation is singular at r = 0. The result is, denot¬ 
ing by r® the a - th gauge group generator, A 4 s A“7’“ 

A 4 = (7) 

8 

B ^r->oo P[ — 7 : + Dirac - S tying] ( 8 ) 

gr^ 

The Abelian projection appropriate to expose stable 
monopoles in QCD is the gauge 54 A 4 = 0 and A 4 diag¬ 
onal. This gauge is perfectly defined and has no gauge 
freedom left except the residual t/(l). 

A 4 in this gauge can be traded with the Polyakov line 
L{r), which is defined as the trace of parallel transport 
along the time axis with periodic boundary conditions 
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in time, is gauge covariant and hence its trace is gauge 
invariant. 


L{r^ = Pexp( r igA 4 (r,T)dT) (9) 

Jo 

Here P denotes ordered product, T is the temperature or 
j the extension of the system in the temporal direction. 
QCD will be defined as the limit for T ^ 0. 

In the gauge 54 A 4 = 0 Eq.Q simplifies to 


,.gA4{r) 

L(r) = exp(i ) 


( 10 ) 


The gauge in which A 4 is time independent and diago¬ 
nal coincides with the gauge in which the Polyakov line 
is. The appropriate abelian projection for monopoles 
as stable classical solutions is the one in which the 
Polyakov line is diagonal. 

Notice that Eq.’s Q are invariant under the change 
A, ^ A, A 4 —> -A 4 , and, as QCD, under parity 
and charge conjugation. By invariance under parity [ 
Ai{r) -Ai{-r) , Ai)(r) — » Ao(-i^] for each solution 
there will be a solution with the same energy having the 
same magnetic field B and opposite electric field -E 
and the same boundary condition at r —» 00 . 

Under charge conjugation, in the usual representation 
for SU(2) with Ti, Tj, real and T 2 immaginary the 
charge conjugate classical solution is minus the com¬ 
plex conjugate, which is of course again a solution with 
the same energy. Changing A 4 ^ -A 4 , which is also a 
symmetry, we get a solution with the original boundary 
conditions, the same electric field and opposite mag- 
netic field. The energy (mass) of the monopole is lllhll 


M = ^u. 

g ^ 


Everything is easily generalized to a generic gauge 
group, but we will only consider S U{N). 

Eor S U{N) there exists one S U{2) subgroup for each of 
the positive roots of the Lie algebra, and with it a 

static monopole solution. For the (N-l) simple roots a, 
, i = - 1 ) the solutions read, in the unitary gauge, 

Af ( 11 ) 

Here T® is the third component of the SU{2) sub¬ 
algebra spanned by the generators corresponding to the 
i - th root , C = 2,C' is the sum of the fundamental 
weights corresponding to the i - th simple roots. The 
space components A® have the same form as in Eq® 
with the index a running in the appropriate S U{2) sub¬ 
algebra corresponding to the i - th simple root. Recall 
that [C',7’^'] = 0, [C\Ti] = = ^'■'[7’', rj,]. This 


makes C -T'^ commute with Tj., which is the condition 
for (fTTT i to be a solution of Eq Q. Invariance under Weil 
transformations fixes the scale yu to be the same for all 
simple roots. 

For the roots a(,^k) = a, + ...dt+k-i, t = 1, ...N - k the 
solution with the same boundary conditions reads 

Af'^^ = yu[C - krf *\l - Hm))] (12) 


In this notation k - 1 corresponds to fundamental 
monopoles. The mass of the monopole (t,k) is k 
times the mass of the monopole for a simple root, 
namely - k^B- The factor k is dictated by the 

requirement that the boundary condition at large r be 
the same as for the fundamental monopoles. Since 
[C, the factor 

k is needed to make [C - kT'^’'‘\ - 0 and have a 

solution. There is again only one free parameter, the 
scale fi. 

Besides the above monopoles there exist their 
excitations ll2ll] . Since yu is arbitrary, it can be changed 
to yu„ = yU -H nAjrT with n an integer and the solution 
Eqdlll to Af ’ = C[yu -H n47TT] - krf V(1 - H(k^„)), 
with = r^n- A gauge transformation 


U)(x4) - exp(i47rn7’C.X4) (13) 

periodic in .X 4 , transforms A^’*' —> A^’^* - CnAn^, or to 
yuk) _ A monopole with the 

same quantum numbers, the same boundary condition at 
large r and mass + ikArm)^. The classical weight 
exp(-5) = exp(-^) is depressed by a factor exp(^^^) 
with respect to original monopole. 

Another series of monopole solutions can be obtained 
by use of the excitations of the charge conjugate solu¬ 
tions. Define p. by the equation 


k/u + {N - k)ij - 2nT (14) 

The charge conjugate of the solution Ea (fT^ is again a 
solution, and, since // is arbitrary, k/j can be changed to 
{N-k)p . The resulting solutionis Cii+{N-k)pT^''^\\- 
H(^{N - k))) again compatible with the boundary con¬ 
ditions at spatial infinity. Indeed it is easily seen that 
C = ^ ^ 0. By use of 

EqdHli Cb + {N- k)T^^’^'>p = yuC'^ -H 2nTTf'‘\ The last 
term can be gauged away by a transformation periodic 
in time of the type Ea (fT3l) . ^ s pr. 

In summary, neglecting excitations, in SU(N) gauge 
theory, there are monopoles (f, k), f = 1, ...N - k 
and k - 1, ..N-l, with scale parameter yuk, independent 
of t, N - k for each value k, to which additional k can 
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be added with scale parameter fik, defined by the rela¬ 
tion fik + iJ.(N - k) — 2nT , for a total of - 1. The 
corresponding A4’s are those of Ea lfTSli for the first set 
t - 1...N - k and, for the second set f = A7 - ^ -H I...N 


Af > \l - H{km (15) 

The generators are A7 x A7 matrices. The are diag¬ 
onal with diagonal elements 
[Tf\,n = ^n,i - S„,,+k for t<{N- k), and 
- 5n,,^k-N for t>{N- k). 

[Tt*\n, = for t<(N-k), 

lT+'‘\m = 5i,6m,t+k-N for t>(N -k). 

It is a general theorem that for a soliton the trace of the 
Polyakov line at large distances is constant on the sphere 
at infinity, independent of the direction lE3l . That value 
is called Holonomy, ^ T r[L(cx))] . In our case it can be 
explicitly computed giving 


1 1 1 sin(^) 

-rr[L(oo)] = -rr[exp(/Ca)].-^ (16) 

Here a = j and 


C = diag{ 


N-l 

2 


N-l 

2 


1 ,... 


A7- 1 
2 


) 


(17) 


The holonomy Ea lfTfil) vanishes for a - ^ and is = 1 
for a - 0 mod ^. 


3. The Polyakov line and Confinement. 

Eor pure gauge theories (no quarks) the trace of the 
Polyakov line averaged over the position (L) is the or¬ 
der parameter for confinement: it vanishes in the con¬ 
fined phase, and is non zero in the deconfined one. 

f 

{L) - exp(--^) where T is the temperature, and Fq the 
free energy of a static quark: (L) = 0 means diverging 
Fq i.e. confinement. 

The symmetry involved is the invariance under the Cen¬ 
tre of the gauge group, which leaves the ground state 
invariant in the confined phase, and is spontaneously 
broken in the deconfined one. In presence of dynami¬ 
cal quarks that symmetry does not exist. Therefore that 
symmetry can not be the symmetry responsible for con- 
finememt, as discussed in Sect[T] since quarks do exist 
in nature. However it can be related to it, and coincide 
with it in the quenched case as we shall argue below. 
Empirically the Polyakov line still works as an order 
parameter: it has a rapid variation at the transition, or a 
peak in the susceptibility, which obeys scaling with the 


appropriate critical indexes as in the quenched case. 
Prom the definition Ea. ifTOl) it could seem that the 
Polyakov line has infrared problems when the tempera¬ 
ture T tends to zero, the limit which defines QCD. In 
fact as r —» 0 the time extension Ajc 4 of the system 
becomes larger and larger, much larger than the corre¬ 
lation length , which is typically of the order of 1/m, 
approximately equal to the inverse of the temperature of 
deconfinement. Since at time differences larger than the 
correlation length the boundary conditions become ir¬ 
relevant, one can conclude that as Ax4 —> oo, (L) a; (L)^, 
k an integer, which implies either (L) = 0 (confinement) 
or (L) = 1, and that, below the deconfining transition, 
Ax 4 = in Eg. ([Tol l can be replaced by If (L) has to 
be the holonomy of monopoles, in Ea. ifTfil) T is replaced 
by A, a = and /j. - ^ A. The size of the monopoles 
is fixed by the correlation length. 

In the confined phase a = ^ in Ea (fT4l l gives p. - fi. 
Monopoles appear in sets of N with the same mass: 
the A - 1 fundamental monopoles degenerate with the 
charge conjugate of the monopole corresponding to the 
root have scale parameter p; The 

N - 2 monopoles corresponding to the roots a' + 

have scale parameter 2p and are degener¬ 
ate with the charge conjugates of the monopoles corre¬ 
sponding to the roots a' +d^+..d^~^ and d^+d^ + ..d^^^, 
and so on till the monopole a' + together with 

the charge conjugates of the fundamental monopoles c?' 

which have scale (N - l)p. The fundamental 
monopoles will dominate since the weight is propor¬ 
tional to ofexp(-^). 

One can compute the Polyakov line distribution of the 
monopoles: the computation is tedious but trivial since 
the exact form of the monopole field is known. The 
contribution of a monopole (t,k) sitting at the origin 
depends on the distance from the centre and for f = 
1,{N - k) has the form 

- Lao = ^exp(-/a[ ^ ^ ^ - f]) 

cyJc cyk 

[cos(—))-cos(—H(kmm 

Eor t - N - k + 1, ...N 

„ , 2 N + I ak 

C’‘‘Hr)-Lao = — exp(-!a[^-fjjexpiy 

cyk cyk — 

[C0S(y)-C0S(yH(k^))] (19) 

In the confined phase a - a - = ^ and the 

two equations (fTSl l and (fT^ have the same form, with 
the index t running from 1 to A, namely 

'5'JT 

C’^HP) - Lao ^ exp(it—)<i>ik, ^) (20) 
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where 

2 N-k+1 

= ^exp(-!a'[- - -]) 

cyk cyk 

[cos{—)-cos(—H(k^))] (21) 

In the confined phase Loo = 0 rEa lfTfil ll and the contribu¬ 
tion of the monopole t, for whatever value of k, is pro¬ 
portional to the t - th N root of 1. If the weight is inde¬ 
pendent on t the contributions of the monopoles sum up 
to zero. It is easily seen that the argument holds for all 
higher excitations of the monopoles. Adding monopoles 
to the confining vacuum leaves (L) locally unchanged. 

In the simple case of SU{2)k - 1 and Ea (l20l i has the 
form - exp(/7rf) cos(|[cosh(f) - |]). The two val¬ 
ues f = 1 and t - 2 give opposite sign: if the distribution 
is the same: this makes the average value (L) = 0. 

4. The vacuum state. 

If the ground state in the confined phase in the classi¬ 
cal limit is a field configuration of monopoles with com¬ 
patible boundary conditions at large distances, the size 
of the monopoles is fixed to be yu = ^A. The shape 
of the monopoles Eg. (1211 1 and their distribution deter¬ 
mines the distribution of the Polyakov line in complex 
space. We shall consider monopoles with k = 1: the 
result is trivially extended to higher mass monopoles. 
The contribution of a single monopole configuration is, 
apart from a factor proportional to the A-th root of the 
identity [Eq.s (l20t and (l2T1 i ], 

L(^) = |[cos(^//(^)) - cos(^)] (22) 

H{^) = coth ^ - I grows monotonically from 0 at ^ = 0 
to 1 at ^ » oo, and L decreases correspondingly from 

Tmax - cos(^) -I- 1] to 0. At Small values of ^ 

H(^) » I and 

(23) 

Since ^ a; -^(^)^, the distribution in space Pi^)d^ - 
can be traded with the distribution in L 
P{L) - L’(^)/|^| . If Nm is the number of monopoles 
in the volume V and pM their density each of them will 
contribute to the distribution near L,„„^ in the same way, 
giving 

P(L) = A^(A[L„,,,-L])5 (24) 

with K - The value a - ^ has been used. 

Eor S U{2) the distribution observed on the lattice is 


P{L) Si |(1 - L^)5 and the number of monopoles in a 
cube of size of the order of the correlation length is of 
order unity. A systematic comparison of Ea (l24li with 
the distributions observed on the lattice at various N 
would be of interest. 

What we have done so far is to take seriously the idea 
that monopoles in QCD exist as stable topological 
configurations and draw the consequences: the first 
consequence of this assumption is that the monopoles 
belong to the abelian projection in which the Polyakov 
line is diagonal. In addition they are BPS monopoles, 
i.e. self-dual or anti-self-dual configurations. This 
feature can in principle be verified in lattice configu¬ 
rations: an excess of magnetic flux trough a plaquette 
(i,j) of an elementary cube should correspond to an 
excess of electric flux through the space-time plaquette 
(0, k) with k - i A j. 

We have argued that, below the deconfining transition 
down to zero temperature the effective time extension 
which determines the Polyakov line is the correlation 
length of the system Requiring the holonomy to be 
that of confining vacuum, namely TrL - 0, fixes for 
S U{N) to Y the ratio ^ of the scale parameter of the 
monopoles to the physical scale A of the system. 

Notice that any gauge transformation periodic in X4 
with period j leaves the trace of the Polyakov line, 
unchanged. 

Requiring the compatibility of the boundary conditions 
for all the solutions fixes the mass spectrum of the 
monopoles. The spectrum is specially simple if the 
holonomy vanishes (Confinement): the monopoles 
group into A-plets and the contribution to the Polyakov 
line of their members differ by a factor A-th root 
of the identity, so that they sum up to zero if their 
distribution is symmetric. Adding a monopole leaves 
locally unchanged the Polyakov line, and in particular 
the confining vacuum. 

This picture nicely fits the arguments of Refi ll^ on the 
universality class of the deconfining transition. 

Eor A > 3 pure gauge the transition is first order, 
and the corresponding temperature Tc ~ A: below Tc 
a - ^ (L) - 0 Ea lfTfil l. ft - p. In the deconfined 

phase /i = 0, A4 = 0 and there are no stable monopoles. 


What is observed in lattice QCD configurations |124 , 


25h is that, in the confined phase, stable monopoles 


do exist in the maximal abelian gauge, which wrap 
the time direction of the lattice via periodic bound¬ 
ary conditions, the so-called ’’long monopoles”. These 
monopoles dominate physics, the short ones are irrele¬ 
vant. 
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It is also known that for static monopoles the unitary 
gauge and the maximal abelian gauge coincide j^ . 
Our stable monopoles could then be identified with the 
long monopoles. 

Magnetic charge is conserved, so that monopoles of any 
kind describe closed loops in space-time, most of which 
close without crossing the lattice in the time direction. 
Some of them could be lattice artefacts, coming from 
the method of detection fioll . Empirically all of them 
seem to be irrelevant to physics. 

A guess for the ground state |G) could be a coherent 
state of stable monopoles 

|G) = exp(J3bl)\0) (25) 

with br - brip - 0) the destruction operator at momen¬ 
tum p = 0 of the r-th monopole with k - irEa.s(fT2li. 
(inii ]. To guarantee invariance the operator acting on 
|0) in Ea (l25T l has to be multiplied by its space reflected 
and the result by its charge conjugate. 

A state of the form Ea (l25T l is left invariant, apart for a 
trivial factor, by applying the sum 'Zrb]-, since the weight 
is simmetric in r, [ Ea (l20l i ]. In principle a linear com¬ 
bination of monopoles is not a monopole, since EqQ is 
non linear; it is not trivial to define a Hilbert space for 
monopoles. However a scalar product can be defined in 
terms of the classical field configurations,, which makes 
monopoles with different charges orthogonal, by use of 
the conserved magnetic current. Physically screening of 
the magnetic field produced by dual superconductivity 
makes monopoles practically non interacting. 

Monopole condensation in this picture is Bose conden¬ 
sation as discussed in Ixh . The order parameter is 

p = {G\bM\G) (26) 


The vacuum Ea (l25T l is only a plausible guess. Proving 
it is not easy. 

An alternative approach starts from Calorons which are 
a generalisation of the usual instantons to the case of a 
non trivial Holonomy ll^ ll^ . In Ref| 281 it was shown 
that they have monopoles [dyons] as constituents. Their 
weight can be computed as for ordinary instantons ll^ . 
This allows to compute the weight of the monopoles and 
to circumvent the infrared problem that would emerge in 
computing it for a single monopole, which is charged. 
A generalisation of the instanton ensemble model to the 
case of non trivial holonomy. Even in that model the 
attempt to show that there is a minimum at (L) = 0 in 
the effective action for the Polyakov loop was not suc¬ 
cessful ll^ . The assumption is that all the monopoles 
can be considered as constituents of calorons. Lattice 
data llslll show that this may not be the case; in S U{3) 


gauge theory with no quarks adding monopoles to field 
configurations changes the number if instantons , one 
for each monopole anti-monopole pair. Moreover it is 
known that instanton ensembles have difficulty in ex¬ 
plaining confinement. 

Dual superconductivity of the vacuum can be the way to 
explain confinement in terms of symmetry, and in par¬ 
ticular an alternative way to look at the centre symmetry 
in the pure gauge case. Stable monopoles seem to play 
an essential role. More work is needed to put it on safe 
grounds. 
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